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Abstract 

We introduce a gaussian probability density for the space-time distribution of 
wormholes, thus taking effectively into account wormhole interaction. Using a mean- 
field approximation for the free energy, we show that giant wormholes are proba- 
bilistically suppressed in a homogenous isotropic "large" universe. 
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1 Introduction 



Some years ago it was observed [Q] that topological fluctuations of space-time provide 
a mechanism which can in principle fix the values of some fundamental constants 
of nature. More precisely, it appeared that taking into account the possibility of 
creation of tiny wormholes connecting distant regions of space-time has the effect of 
modifying the coupling constants contained in the model and to provide a probability 
distribution for them. In particular, it was shown by Coleman that, under certain 
assumptions, the probability distribution for the cosmological constant was infinitely 
peaked around zero, thus providing a candidate solution, at least neglecting one-loop 
effects [0, to the long-standing cosmological constant problem. 

The main characteristic of wormholes is that the probability of their creation 
does not depend on the positions of the endpoints xq and x\, at least when the 
two points are distant. In a first approximation, one can think that the wormholes 
are dilute, so that \x\ — Xq\ is in the average great with respect to the wormhole 
size a, and one can forget about the wormhole interaction which is likely to appear 
when the endpoints become close to each other. This approximation was assumed 
in all previous analysis of this subject. However, it was evident from the very 
beginning [|J that this simplified model ran over some inconsistencies, like the so- 
called giant wormhole paradox, that is, the non-physical prediction that wormholes 
of macroscopic sizes should be generated. 

In this paper we introduce a finer description of the process of wormhole creation 
and destruction, and show that in this case giant wormholes are probabilistically dis- 
favoured. More precisely, we assume that the wormhole space-time distribution be 
described by the gaussian probability density Q2.3|) , which corresponds to a suppres- 
sion of the probability of creation of a wormhole when the two endpoints come close 
to each other. We then evaluate the functional integral for the theory using methods 
which were previously introduced by one of us [|J to deal with similar problems in 
solid state physics. Using the mean-field approximation (|3.12 ) for the free energy, 
we show that the free energy for the theory is proportional to A 4 a 2 /£p (where A is an 
energy scale < 1/fp and a is the effective size of the wormholes), and therefore has a 
minimum for small values of a. This is our main physical result, which implies that 
giant wormholes are probabilistically suppressed in a homogenous isotropic "large" 
universe. This phenomenon is a large-scale, low-energy "cooperative effect" between 
the wormhole density ( [2.1|) and space-time distribution ( |2.3|) . Previous analysis of 
the problem neglected a self-consistent determination of the wormhole size a because 
they did not take into account the wormhole space-time distribution, arriving this 
way to the cited giant wormhole paradox. Furthermore, our analysis shows that 
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the result is quite independent from the "microscopic details" of the gravity-matter 
interaction, at least in the domain of validity (1/a < A) of our model. 

2 Definition of the model 

As a starting point, let us assume that we have performed integration over short- 
range topology fluctuations. It is natural to suppose that a reasonable trial wave 
function, after eliminating such degrees of freedom, will give us some effective theory 
of geometrical and matter degrees of freedom on the wormholes background. More- 
over, it seems reasonable that the only trial parameter left in the theory should be 
the effective size of the wormholes. Thus, let us assume that the probability density 
for the creation of a wormhole of size a be of the form 

w = e~ ca2 (2.1) 

where c = const = 0(1) and the length a is measured in units of the Plank length 
ip. 

Now, it is reasonable to think that the above wormhole density distribution for 
finite values of a affects also the space-time distribution of the wormhole itself. 

Of course, it is difficult to build a complete theory of the microscopic dynamics 
of wormholes, but it seems reasonable as a first step to describe this dynamics self- 
consistently by taking into account the size a both in the density and in the space- 
time correlations. Let us in the following formulate a simple model of universe with 
wormholes in which one may calculate explicitly the effect of the size a, and the 
main result shall be that giant wormholes are probabilistically suppressed. 

For this purpose, we describe the space-time distribution of the "ends" of the 
wormholes by the local field 

a(x) = aS(x — x Q ) (2.2) 
and we postulate the following weight function for the field a(x): 

V[a(x)\ = exp(— - J dxdx a(x)a(x)K(x — x')), (2.3) 

Here the correlation function K(x), which is the amplitude for a wormhole insertion, 
is assumed to be gaussian: 

K(x) = K exp(-ax 2 ), (2.4) 

where K ~ exp(— S w ) and S w is the (constant) wormhole action. The ansatz ( |2.3| ) 
allows us to take effectively into account the finite size of wormholes in the space- 
time distribution. As a matter of fact, if we substitute in (|2.3| ) the distribution of 
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two point-like wormholes 



a(x) = ao5(x — Xq) + a±S(x — x\) (2.5) 



we see that the distribution ( |2.3| ) describes a suppression of probability when \x\ — 
Xq\ < a, and we get factorization when \xq — x\\ — > oo. 

As a first order of approximation one may think that the field a, which describes 
the large-scale structure of the space-time distribution of wormholes, is coupled to 
gravity by volume effects only. In other words, we may assume that the microscopic 
structure of the gravitational field is not relevant to describe the size distribution of 
the wormhole dynamics. Therefore the gravitational degrees of freedom associated 
with the spin-2 part, i.e. h^ u = g^ v — (g^), which control the small-scale structure, 
can be neglected at this level. As a consequence of this picture, we impose the 
following coupling of a(x) with gravity: 

Sint = J d A x\Jg{x) a{x), (2.6) 

where the "volume field" yjg is expanded as 



yfgiti = l + tp{x), M<1. (2.7) 

The model that we have in mind is given by the path integral 

Jva(x) e W {-\fd^x'a{x)K{x-x')a{x') + sj) 

= - 1 exp f~ f d 4 xd 4 x' ip(x)K~\x - x')ip(x')) (2.8) 
V det K V 2 J ) 

where ( p.8| ) is defined up to a normalization and we have excluded the term linear 
in (p. Indeed this contribution may be neglected if we integrate ( |2.8j ) over </?, since it 
represents an arbitrary fluctuation field around the vacuum ^y 7 ^) = 1. Of course, 
looking in this way, we are taking into account quantum gravity volume effects with 
a functional measure given by exp(— S(ip)), where S((p) is the effective action 

S((f) = - | dxdx (p(x)K~ 1 (x — x')ip(x'), (2.9) 

coming from the above functional integration over a. We get Coleman's picture 
when we consider only the global degrees of freedom of the field a(x). For instance, 
if we substitute into the Da-integral the factor 

5 (a - J d 4 x a(x)) (2.10) 

and perform the Da-integration, we would have a theory with the global variable a 
playing the role of an effective wormhole collective coordinate in Coleman's spirit. 
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In any case, the effective model based on ( |2.9| ) is only ultralocal for <f. In a 
realistic scenario of quantum gravity coupled with matter fields we must investigate 
also the effect of the density and space-time wormhole distribution over the coupled 
system matter-gravity. So, a less crude model is obtained if we add to the non-local 
action ( |2.9| ) a matter field contribution S mat represented by a massless Klein-Gordon 
field coupled to the geometry g^ u . Assuming again that only volume effects are 
relevant, we may replace the covariant Laplacian in the S mat with the flat one, and 
using |T7l ), S mat looks as 



5 ma t = / d 4 xp(x)(l + <p(x))Ap(x) (2.11) 



where A is the flat, Euclidean 4-dimensional Laplacian. 

The coupled model is now described by S(<f) + S mat (p, ip). Performing firstly the 
p-functional integration, we get a free energy of the form (in Euclidean signature) 

F[^]~ilogdet[(l + ^(x))A]. (2.12) 



Of course we must understand the functional determinant in fl2.12 ) as a renormalized 



one. For this aim in the following section we shall adopt a Schwinger-DeWitt proper- 
time scheme, together with an "-^ expansion". 

3 Mean Field Approximation 

The complete quantum partition function of our model is given by 

Z = JVf exp (-- J d^xd^x' f{x)K-\x - x')f(x') - F[<p]\ 

= {exp{-F[p])) v (3.1) 

where the average symbol means 

(0) v = JVfO exp J fK~ l ^J (3.2) 

At this point we use the mean-field approximation 

Z ~ exp (- (F[<p]) v ) (3.3) 

Let us remind that, due to the convexity of the exponential function, one has the 
Peierls inequality 

(e A ) > e< A >. (3.4) 

Here we assume to take the lower bound of ^e" 4 ^), since in our case A is given by a 
functional — F[ip] of tp, and the kernel ( |2.4| ), goes very rapidly to zero for \x— y\ — > oo; 



this is just the large-range approximation underlying all our model. 
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In other words, the true free energy T = — log(Z) is approximated by 

T * (F[cp]) v = (I logdet[(l + <p)A]) . (3.5) 
By using now the Schwinger-De Witt proper-time formalism, we may write 

logdet[(l +<p) A] 
= - — Tr exp (— t[(l +(p)A]), e -> 0+ (3.6) 



where Tr stands for the space time trace of the heat-kernel operator in ( |3.6| ) and 
can be represented as a Schrodinger-like path integral over closed paths according 
to 

Tr exp (-*[(1 + if) A]) (3.7) 
= [d 4 x [ Vx(r) exp (- ( dr-(l + (p[x(r)}) x{t) 2 ) . 

J Jx(0)=x(t)=x V JO 4 / 

Putting ( |3.6| ) and (3.7) in (|3.5| ), we get 

T ~ lim + y U(pexp(^—-J dxdx'<p(x)K(x — x')(p(x') S j ■ (3.8) 



d 4 x I — [ Vx(r)expl — f f-(l + (p[x})x 



t 



t Jx{0)=x(t)=x { JO \4 

At this point, the functional integration can be done j3J, giving the expression 

T ~ lim -- / d 4 x [°° - I Vx(t) ■ (3.9) 

e^0+ 2 J Je t Jx(0)=x(t)=x 

■ exp [ - [ t dr^^ + - f f drdr' x(t) 2 K[x(t) - x(t')]x(t') 2 ) 
\ Jo 4 2 jo Jo / 

It is important that the path integral over T>x{j) is translationally invariant and 
that the integration over d 4 x gives simply the volume factor V (actually a covariant 
volume with respect to a background metric such that yfg = 1): 

/ = lim / — / Dj(t) ■ (3.10) 

e->0+Je t Jx(0)=x(t)=0 

■exp I- f ' dr^- + \ f* C drdT'x(r) 2 K(x(T) - x(t')x(t') 2 )\ . 
[ Jo 4 4 Jo jo J 

Path integrals of this kind can be found in the theory of disordered systems. 
Unfortunately, it is impossible to calculate exactly even this very simplified path 
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integral. However, we can estimate it in some sense non-perturbatively and obtain 
rough numerical estimates and an enough good qualitative picture. 

Really, our x(t) lies in 4-dimensional space, here assumed isotropic and homoge- 
nous, and it seems to be reasonable that the main contribution to the path integral 
( |3.11[ ) is given by trajectories having high degree of isotropy. For such trajectories: 



< x(t)x(t') > = 0, 
< (x(r)x(r')) 2 > = — < x 2 (r) > 2 (3.11) 

where D is the dimension of the space. In our case D — 4, and we can formulate 



the mean-field approximation for the path integral ( 3.11 ) by requiring that: 



K(x(t) — x(t')) = ca 2 exp^ — -(x(r) — x(t')) 2 ^ 



« ca 2 exp{-\x 2 (T))exp(-\x 2 (r')) . (3.12) 



a 2 / V a 2 



(It is important to notice that this decomposition becomes exact in the limit a — > oo; 
see below). In this approximation the calculation of Q3.ll ) can be reduced to solving 
the quantum mechanical problem 



• / Vx(t) exp < — dr 

Jx(0)=x(t)=0 ' Jo 



, , CO rjf f + OO „ „ 

lim / — / dse~ s /ca G 8 {t\0,0), (3.13) 



oo 



+ 0+ Je t 

where G s (t\x,x') is the Green function 

d t G s (t\x,x') = ^(l + 4se- x2/a2 ) AG s (t\x,x'), (3.14) 
G s {0,x,x r ) = 5 {4) {x-x'). 

Some remarks are in order. In our mean-field approximation the fluctuations of 
gravity are described by the global variable s, which has the meaning of a sort of 
order parameter. Indeed the parameter s can be considered as an amplitude of the 
volume-field ^fg fluctuations. It seems that the very existence of an order parameter 
of this kind has more general meaning and more deep origin than our approximation. 
We see from ( |3.15|) that at s = —1/4 we will have G s (t\x, x') very singular at x — > x'. 
Such fluctuations can be described as bags in space-time, where ^fg — > 0. 

On the other hand, the domain s < —1/4 corresponds to fluctuations with 
^fg < 0, i.e., we have a first-order phase transition, since also the covariant volume 
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factor V = / d 4 x \fg in ( |3.11| ) changes sign, and hence we pass from a regime, say, 



with JF > to one entropy dominated, i.e. T < 0. From a quantum-mechanical 
point of view, a negative free energy should imply vacuum decay, so that in order 
to avoid this possibility we restrict ourselves to the domain of integration s > —1/4 
(i.e., ^/g> 0). 

We come thus to the final formulation of our model: the volume density of free 
energy has the form: 



sx 2 e x ' 2 / a2 



Jx{0)=x(t)=0 { Jo \ 4 

In (|3.16l ) is understood the usual zero-point energy subtraction of Tr log A = f . 



4 General Estimates 



A limiting case in which (|3.16| ) can be computed is obtained by setting a — > +oo. 



In this limit we can neglect the x-dependence of the correlator, i.e. we can put: 



e~ x2/a2 ~ 1 (4.1) 



Then, one gets 

/- lim [ + ™- x p^[ +OO dse- s2 / ca2 -<0\e- t{1+ ^ A \0>-f (4.2) 

e^0+ J e £ V KCa 2 7-1/4 

This quantity requires an ultraviolet regularization. Let us use a momentum cut-off 
regularization by defining 

< 0|e-'( 1+4 *) A |0 >= ^4 [ A k 3 dke- tk2 ( 1+4s \ (4.3) 
(27r) 4 7o 



where A is the UV-cutoff. By (|4.3|) , the leading contribution to ( |4.2| ) in A is approx- 
imately equal to (see Appendix): 

/ « -AW8(1 - « - 2 A 4 ca 2 . (4.4) 

Since in our model we have neglected systematically all quantum gravity short-range 
effects, we may assume that the energy-scale A is upper-bounded by the mass Planck 
scale l/£p (in c = % = 1). 
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5 The General Case 



For arbitrary a the quantum mechanical problem ( |3.16| ) cannot be solved exactly, but 



we see from the previous computation that the main contribution is given by high- 
momentum fluctuations and that we can use semiclassical formulas. It is convenient 
to use Laplace (or Fourier) transform. 

The Laplace transform of the Green function (|3.15| ) in the t-variable has a cut 



along the imaginary semiaxis; the inverse Laplace transform can be reduced to an 
integral of the jump of G on this cut over this semiaxis. The resulting expression, 
after (it-integration, has the form: 



'-if 

where G s (u + iO\x, x') obeys the equation 



/-/ —logo;/ dse- s ' ca \ -.ImG s (a; + iO|0,0)-/o, (5.1) 

jo 7T 7-i/4 V nca 2 



(V(x)A + u)G s (u\x, x') x — x'); (5.2) 

V(x) = 1 + tp(x) 

The Green function G s (u\x, 0) is spherically symmetric: 

G 8 (u\x, 0) = G 8 {u,r), (5.3) 

and taking into account the small imaginary part of to we obtain its semiclassical 
limit in the form: 

G.{u,r) ~ ([dr'k{r')) , (5.4) 

Here we have set 



and Hi is the Hankel function of first kind. It is easy to see that the value of 
Im (? s (u;,r — > 0) is the same as in the case a —>■ oo and that we obtain the same 
result (|4.4|) . As a consequence the main contribution would be given by the domain 
k(0) < A, i.e. st > 1/A in the notations of the Appendix. It should be stressed 
that the mean field approximation becomes exact in the limit a — > oo. On the other 
hand, a is large for the dominating fluctuations. Thus, we can hope that the main 
contribution to the free energy in the giant wormhole limit is taken into account by 

(B. 



The cutoff A has physical meaning and is defined by quantum gravity effects. As 
we have noticed above it is natural to think that A ~ l/£p- Our estimations have 
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relation to the case that l/a < A < 1/tp only. In this limit we see that the volume 
density of free energy increases with a as A 4 a 2 , and the energy minimum is realized 
at small a. 
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6 Appendix 

We estimate the integral (O), ([y|. Let s = \{s - 1), then: 



f * iW r 7 r d§ exp (- ieb (5 - 1 



l_ e -¥A 2 )_^A 2 e-T A n -/„. (6.1 



sH 2 iv / 4 

We divide the domain of dt (is- integration in two domains, st < 1/A 2 and st > 1/A 2 . 
The first domain gives a contribution which cancels the fo, while the second gives 



1/2 /■+« eft r+°° / 1 2 \ 1 

1/A 2 <ft / l / 2 \ 4 2 

l/(4A2vW) ^ V 4A2t V 7TC0 5 J ~ AW' > ' ^' 2 ' 

The domains of (it-integration 1/A 2 < t and t < l/(4A 2 v / co 2 ) give a small contribu- 
tion. 



10 



References 



[1] S. Coleman, Nucl. Phys. B307 (1988) 867; 

S. Giddings and A. Strominger, Nucl. Phys. B307 (1988) 854. 

[2] J. Polchinski, Phys. Lett. B 219 251; 

A. Carlini and M. Martellini, Phys. Lett. B 276 36. 

[3] I. Klebanov, L. Susskind and T. Banks, Nucl. Phis. B317 665; 
W. Fishier and L. Susskind, Phys. Lett. B 217 (89) 48. 

[4] I. Kolokolov, Phys. Lett. 123 (87) 302. 



11 



